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The periodic response of a metallic or a superconducting ring to an external magnetic flux is 
one of the most evident manifestations of quantum mechanics. It is generally understood that 
the oscillation period hc/2e in the superconducting state is half the period hc/e in the metallic 
state, because the supercurrent is carried by Cooper pairs with a charge 2e. On the basis of the 
Bardeen-Cooper-Schrieffer theory we discuss, in which cases this simple interpretation is valid and 
when a more careful analysis is needed. In fact, the knowledge of the oscillation period of the 
current in the ring provides information on the electron interactions. In particular, we analyze the 
crossover from the hc/e periodic normal current to the hc/2e periodic supercurrent upon turning 
on a pairing interaction in a metal ring. Further, we elaborate on the periodicity crossover when 
cooling a metallic loop through the superconducting transition temperature T c . 
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I. INTRODUCTION 

One of the most important properties of superconduc- 
tors is their perfectly diamagnetic response to an exter- 
nal magnetic field, the Meissner-Ochsenfeld effect. It is a 
pure quantum effect and therefore reveals the existence of 
a macroscopic quantum state with a pair condensate. A 
special manifestation of the diamagnetic response is ob- 
served for superconducting rings threaded by a magnetic 
flux: flux quantization and a periodic current response. 

Persistent currents and periodic flux dependence are 
also known in normal metal rings and best known in 
form of the Aharonov-Bohm effect predicted theoreti- 
cally in 1959^. Since the wavefunction of an electron 
moving on a ring must be single valued, the phase of 
the wave function acquired upon moving once around 
the ring is a integer multiple of 2ir. A magnetic flux 
threading the ring generates an additional phase differ- 
ence 2np = e/(hc) § c &v ■ A(r) = [2ire/(hc)\ ip, where C 
is a closed path around the ring and A(r) the vector po- 
tential generating the magnetic flux tp threading the ring. 
Here, e is the electron charge, c the velocity of light, and 
h is Planck's constant. Thus, the electron wave function 
is identical whenever (p has an integer value and there- 
fore the system is periodic in the magnetic flux <p with a 
periodicity of 



the flux quantum in a normal metal ring. In particular, 
the persistent current J(ip) induced by the magnetic flux 
is zero whenever $ = ip/^o is an integer. 

The periodic response of a superconducting ring to 
a magnetic flux is of similar origin as in a normal 
metal ring, though the phase winding of the conden- 
sate wavefunction has to be reconsidered. On account 



of the macroscopic phase coherence of the condensate, 
flux oscillations must be more stable in superconductors, 
and London predicted their existence in superconduct- 
ing loops already ten years before the work of Aharonov 
and Bohm 2 . London expected that the magnetic flux 
threading a loop is quantized in multiples of $o because 
the interior of an ideal superconductor was known to be 
current free. Although London the pairing theory of su- 
perconductivity was not known yet, he anticipated the 
existence of electron pairs carrying the supercurrent and 
speculated that the flux quantum in a superconductor 
might be $o/2. This point of view became generally 
accepted with the publication of the 'Theory of Super- 
conductivity' by Bardeen, Cooper, and Schrieffer (BCS) 
in 1957-. Direct measurements of magnetic flux quanta 
$o/2 trapped in superconducting rings followed in 1961 
by Doll and NabaueiP and by Deaver and FairbanlP, cor- 
roborated later by the detection of flux lines of &o/2 in 
the mixed state of type II superconductors^. 

It is tempting to explain the $o/2 flux periodicity of 
superconducting loops simply by the charge 2e of Cooper 
pairs carrying the supercurrent, but pairing of electrons 
alone is not sufficient for the $0/2 periodicity. The 
Cooper-pair wavefunction extends over the whole loop, 
as does the single-electron wavefunction, and it is not 
obvious whether the electrons forming the Cooper pair 
are tightly bound or circulate around the ring separately. 
A microscopic model on the basis of the BCS theory is 
therefore indispensable for the description of the flux pe- 
riodicity of a superconducting ring. In this chapter, we 
analyze this problem in detail and focus on a previously 
neglected aspect: how do the <&o periodic flux oscillations 
in a normal metal ring transform into the $0/2 periodic 
oscillations in a superconducting ring? 

A theoretical description of the origin of the half- 
integer flux quanta was first found independently in 1961 
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FIG. 1: Scheme of the pairing of angular-momentum eigenstates in a one dimensional metal loop for (a) $ = and (b) 
3? = <&o/2, as used by Schrieffer irP to illustrate the origin of the $o/2 periodicity in superconductors. Paired are always 
electrons with equal energies, leading to center-of-mass angular momenta q = in (a) and q = 1 in (b). 



by Byers and Yang 9 , by Onsager^, and by Brenig^ on 
the basis of BCS theory. They realized that there are two 
distinct classes of superconducting wavefunctions that 
are not related by a gauge transformation. An intu- 
itive picture illustrating these two types can be found in 
Schrieffer's book on superconductivity, using the energy 
spectrum of a one-dimensional metal ring. The first class 
of superconducting wavefunctions, which London had in 
mind in his considerations about flux quantization, is re- 
lated to pairing of electrons with angular momenta hk 
and — hk, which have equal energies in a metal loop with- 
out magnetic flux, as schematically shown in figure [I] (a). 
The Cooper pairs in this state have a center-of-mass an- 
gular momentum (pair momentum) hq — 0. The pairing 
wavefunctions of the superconducting state for all flux 
values $, which are integer multiples of $o and corre- 
spond to even pair momenta hq — 2fi$/<I>o, are related 
to the wavefunction for $ = by a gauge transformation. 
For a flux value $o/2, pairing occurs between the electron 
states with angular momenta hk and h{—k + 1), which 
have equal energies in this case [figure [l] (b)]. This leads 
to pairs with momentum hq = h. The corresponding 
pairing wavefunction is again related by a gauge trans- 
formation to those for flux values $ which are half-integer 
multiples of <f>o and correspond to the odd pair momenta 
hq = 2M>/$ . 

For the system to be <&o/2 periodic, it is required that 
the free energies of the two types of pairing states are 
equal. Byers and Yang, Onsager as well as Brenig showed 
that this is in fact the case in the thermodynamic limit. 
The free energy consists then of a series of parabolae 
with minima at integer multiples of $0 (corresponding 
to even pair momenta) and half integer multiples of $0 
(corresponding to odd pair momenta). If the arm of the 
ring is wider than the penetration depth A, the flux is 
quantized and the groundstate is given by the minimum 
closest to the value of the external flux. However, in mi- 
croscopic finite systems this degeneracy of the even and 



odd q minima is lifted, although their position is fixed 
by gauge invariance to multiples of $o/2. The restora- 
tion of the $o/2 periodicity in the limit of large rings 
was studied only much latei^ED. We study the revival 



of the $o/2 periodicity in sections II A and II B for a one- 
dimensional ring at zero temperature and investigate the 
effects of many channels and finite temperatures in sec- 
tion nn 

From the flux periodicity of the free energy, the same 
flux periodicity can be derived for all other thermody- 
namic quantities^. A clear and unambiguous observa- 
tion of flux oscillations is possible in the flux dependence 
of the critical temperature T c of small superconducting 
cylinders. Such experiments have been performed first by 
Little and Parks in 1962^3®!. They measured the resis- 
tance R of the cylinder at a fixed temperature T within 
the finite width of the superconducting transition and de- 
duced the oscillation period of T c from the variation of R. 
These experiments confirmed the $o/2 periodicity in con- 
ventional superconductors very accurately. At this stage 
the question of the flux periodicity in superconductors 
seemed to be settled and understood. The interest then 
shifted to the amplitude of the supercurrent and also the 
normal persistent current and their dependence on the 
ring size, the temperature, and disordeiP^BS. However, 
the influence of finite system sizes on the flux periodicity 
remained unaddressed. Earlier, certain experiments had 
already indicated some unexpected complications. E.g., 
Little and Parks pointed out in reference^ that in tanta- 
lum cylinders they could not detect any flux oscillations 
in R at all. Even more peculiar were the oscillations ob- 
served in an indium cylinder where signs of an additional 
$0/8 periodicity were clearly visible^. This was surpris- 
ing because indium is a perfectly conventional supercon- 
ductor otherwise. These results remained unexplained 
and drew attention only years later, when flux oscilla- 
tions of unconventional superconductors were studied. 

In the meantime a new type of flux sensitive systems 
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FIG. 2: (a) Fourier transform J C (T) of the critical current J C (H) measured by Schneider et al. on a 24° grain boundary SQUID 
at T = 77 K as a function of the applied magnetic field where $o/2 = 6.7u.T— (b) Critical current J C (H) over a 24° grain 
boundary SQUID at T = 4.2 K, where $o/2 = 2.7 U.T. Clearly visible is the abrupt change of periodicity at noH w ±5 



was advanced: superconducting quantum-interference 
devices (SQUIDs). The measurement of flux oscilla- 
tions in SQUIDs is similar to the Little-Parks experi- 
ment. Here the flux dependence of the critical current 
J c through a superconducting loop including one or two 
Josephson junctions is measured. This has the advan- 
tage that flux oscillations can be observed at any tem- 
perature T < T c , and they are most clearly visible in 
the critical current J c . SQUIDs fabricated from conven- 
tional superconductors have been used in experiments 
and applications for five decades, and they proved to 
oscillate perfectly with the expected flux period $o/2- 
It was therefore a surprise that flux oscillations with 
different periodicities were found in 2003 by Lindstrom 
et alW± and Schneider et aZP12SI m SQUIDs fabricated 
from films of the high-T c superconductor YBa 2 Cu30 !/ 
(YBCO) where the Josephson junctions arise from grain 
boundaries. Flux trapping experiments in loops showed 
that flux quantization in the cuprate class of high-T c su- 
perconductors occurs in units of <J>o/2-^, identically to 
what has been observed with conventional superconduc- 
tors. In addition, Schneider et al. observed a variety 
of oscillation periods, depending on the geometry of the 
SQUID loop, the grain-boundary angle, the temperature, 
and the magnetic-field range of the SQUID. 

Two distinct patterns of unconventional oscillations in 
YBCO SQUIDs have to be discerned. The first kind con- 
sists of oscillations which have a basic period of $o/2, 
overlaid by other periodicities, such that the Fourier 
transform J C (T) of J c ($) contains peaks appear which 
do not correspond to the period <&o An example for 
such a measurement are shown in figure [2] (a). The peaks 
at integer values of T correspond to higher harmonics of 
$o/2, and their appearance is natural. However, there 
are clear peaks at T = 1/2 (red arrow) and T = 5/2, 



which correspond to $0 periodicity and higher harmon- 
ics thereof. The origin of the $0 periodicity in those 
experiments is so far not conclusively explained. There 
was, however, extensive research on the flux periodicity of 
unconventional (mostly d-wave) superconductors, which 
revealed that the periodicity of the normal state persists 
in the superconducting state if the energy gap symmetry 
allows for nodal states^ESI. This effect derives directly 
from the analysis in this book chapter and is discussed 
in detail in reference^. 

The second kind of unconventional oscillations is more 
intriguing. In several different YBCO SQUIDs, the pe- 
riodicity of sinusoidal oscillations changes abruptly with 
increasing magnetic flux. In the measurement shown in 
figure[2](b), the period is 4>o/4 for small flux, and changes 
to $o/2 at a critical flux. As a possible explanation for 
the appearance of $o/4 periodicity, an unusually pro- 
nounced second harmonic in the critical current J c of 
transparent Josephson junctions was proposed or, more 
fundamentally, an effect of interactions between Cooper 
pairs, leading to the formation of electron quartets^. 
The observation of similar abrupt changes to other frac- 
tional periodicities like $0/6 and $0/8 render this find- 
ing even more striking since it could indicate a transition 
into a new, non-BCS type of superconductivity. This 
concept, which we sketch briefly in the Conclusions, is a 
complex and promising topic for future research on un- 
conventional superconductors. 



II. THE PERIODICITY CROSSOVER 

In this section we introduce the periodicity crossover 
and consider first the simplest model containing the rel- 
evant physics: a one dimensional ring consisting of N 
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FIG. 3: The simplest description of the many-particle state in 
a flux threaded loop, we use a tight-binding model on a dis- 
crete, one-dimensional ring with TV lattice sites, lattice con- 
stant a and radius R = Na/2n. The magnetic flux $ is 
confined to the interior of the ring and does not touch the 
ring itself. 

lattice sites and a lattice constant a (figure [3]). The ring 
is threaded by a magnetic flux $ which does not touch 
the ring itself. We use a tight-binding description with 
nearest-neighbor hopping parameter t, which sets the en- 
ergy scale of the system. We start from the flux periodic- 
ity of the normal metal state of the ring, which varies for 
different numbers of electrons in the ring. On this basis 
we introduce a superconducting pairing interaction and 
investigate the flux periodicity of the groundstate upon 
increasing the interaction strength. For a ring with a 
finite width (an annulus) we investigate the flux depen- 
dence of the self-consistently calculated superconducting 
order parameter and study the temperature driven peri- 
odicity crossover when cooling the ring through the tran- 
sition temperature T c . 

A. Normal state 

The tight-binding Hamiltonian for an electronic system 
including a magnetic field is straightforwardly formulated 
using the annihilation and creation operators Cj s and c| s 
for an electron with spin s on the lattice site i: 

n Q = -tJ2 e^4c is -p^cl*.. (2) 

Here denotes all nearest-neighbor pairs i and j, 

s =t, I- The magnetic field B = V x A enters into 
the Hamiltonian ^ through the Peierls phase factor 

ifij = (e/hc) f? dl ■ A. The chemical potential fi controls 
the number of electrons in the ring. The flux periodicity 
is easiest to discuss for a particle-hole symmetric situa- 
tion with fi = 0, for which the Fermi energy is E-p = 0. 
We will later address the changes introduced through an 
arbitrary fj,. 

We assume that the TV lattice sites are equally spaced 
along a ring with circumference 2irR = TV a (figure uty . 
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FIG. 4: Energy spectrum of a discrete one-dimensional ring 
with TV lattice sites, fj, = and (a) TV/4 is an integer, (b) TV/4 
a half integer, and (c) TV an odd number. In (a) and (b), levels 
cross Ef — at integer (a) or half-integer (b) values of <p. For 
odd AT, two different spectra are possible [TV = 4n + 1 (left) 
and TV = 4n — 1 (right)], and for both, two levels cross Ep 
within one flux period (red points), en denotes the maximum 
value of the Doppler shift. 

It follows that the Peierls phase factor for a magnetic 
field focused through the center of the ring simplifies to 
ifiij = 2irtp/N, where <p = <&/<l>o is the dimensionless 
magnetic flux. The Hamiltonian |2]) is then written in 
momentum space as: 

H = ^2e k (<p)cl s c ks (3) 

k,s 

where ct creates an electron with angular momentum 
hk. The energy dispersion is 

e fc (^) = -2tcosfc^ -fi. (4) 

The eigenenergies depend on the flux only in the com- 
bination k ~ (p, as is shown in figure [4] for three different 
cases: (a) TV/4 is an integer, (b) TV/4 is a half integer, 
and (c) TV is an odd number. The ip dependent shift 
in Ck(f) is known as the Doppler shift since it is propor- 
tional to the velocity of the corresponding electron. In all 
three cases, the spectrum has obviously the periodicity 
1 with respect to tp. However, the flux values, for which 
an energy level crosses are different. This number 
dependence, sometimes referred to as the "parity effect" 
, is characteristic for discrete systems and not restricted 
to one dimension. It was discussed in detail in the con- 
text of the persistent current in metallic loops^BESl and 
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FIG. 5: (a) Total energy E(ip) as a function of the magnetic flux tp. If the number of electrons on the ring N is a multiple of 
four, then E(ip) has minima at integer values of tp (turquoise). If N/4 is an integer, then the minima are at half-integer values 
of <p (blue). If N/4 is a half-integer, the parabolae are shifted by 1/2 (turquoise). The gray lines above the crossing points of 
the parabolae correspond to possible excited states, (b) Persistent current J(ip) corresponding to the systems described in (a). 
The purple curve shows the current obtained for odd N. 



also in metallic nano clusters^ it is also essential for the 
discussion of superconducting rings. 

Physical quantities of the normal metal ring can be ex- 
pressed through the thermal average n s (k) of the num- 
ber of electrons with angular momentum k and spin s: 
n s (k) — (c\. s Cks) = /(efc(v)); with the Fermi distribution 
function /(e) = 1/(1 + e € ^ ksT ) for the temperature T. 
The groundstate is given by the minimum of the total 
energy E of the system 

E(<p) = (H ) = ^2e k (cp)n s (k), (5) 

which is a piecewise quadratic function of the magnetic 
flux. The momentum distribution function n s {k) also 
depends on the magnetic flux only in the combination 
k — ip. The sum over k in equation ([5| directly renders 
the $0 flux periodicity of E(ip). However, the position 
of the minima of E(ip) depends on the highest occupied 
energy level and therefore also shows a parity effect (see 
figure [5]) . 

The energy E(ip) is maximal for those values of ip. 
where an energy level reaches Ep (red points) and has 
minima in between. If N / 4 is an integer, then the min- 
ima of E(<p) are at half-integer values of ip (figure [5] (a), 
light blue curve), whereas if N/4 is a half integer, the 
minima are at integer values of <p (figure [5] (a), dark blue 
curve). If N is odd, two different (but physically equiv- 
alent) spectra for N = An ± 1 are possible, and for both, 
two levels cross E-p in one flux period. This results in 
a superposition of the two previous cases and there are 
minima of E{ip) for both integer and half-integer values 
of ip; E(ip) is therefore <& /2 periodic. 

The normal persistent current J {if) = 
— (e/h)dE(ip)/dip (see equation [14] below) jumps 
whenever an energy level crosses Ep, because the 
population of left and right circulating states changes 
abruptly [figure [5] (b)]. The occupied state closest to 



Ep contributes dominantly to the current, because all 
other contributions tend to almost cancel in pairs. The 
Doppler shift decreases with the ring radius like 1/R 
[c.f. equation E|] and so does the persistent current. 



B. Superconducting state: Emergence of a new 
periodicity 

The theory of flux threaded superconducting loops was 
first derived by Byers and Yangp, BrenigSl, and On- 
sagei^ on the basis of the BCS theory. They showed 
the thermodynamic equivalence of the two superconduct- 
ing states discussed above in the thermodynamic limit. 
However, in a strict thermodynamic limit the persistent 
(super-) current vanishes, and therefore a more precise 
statement is necessary with respect to the $o/2 period- 
icity of the supercurrent. Here we analyze the crossover 
from $0 periodicity in the normal metal loop to the $0 /2 
periodicity in the superconducting loop upon turning on 
the pairing interaction. The discussion of this crossover 
enables precise statements about the periodicity. 

For a one-dimensional superconducting loop (or any 
loop thinner than the penetration depth A), finite cur- 
rents flow throughout the superconductor. The mag- 
netic flux is consequently not quantized, only the fluxoid 
$' = $ + (A/c) § dr • J(r) is, which was introduced by F. 
LondorP. The flux <!> is the total flux threading the loop, 
including the current induced flux, and A = 4ir\ 2 /c 2 . In 
the absence of flux quantization, (p = (e/hc) ip is a con- 
tinuous variable also in a superconducting system with a 
characteristic periodicity in ip. 

In this section we focus on the emergence of a new pe- 
riodicity when a superconducting order parameter arises. 
We therefore include an attractive on-site interaction of 
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FIG. 6: Energy E(<p) in the superconducting state with q = (a, c) and q = 1 (b, d) for N = 50. The upper panels (a, b) show 
the "small gap" case with A = 0.05 1 and the lower panels (c, d) the "large gap" case with A = 0.2 1. 



the general forrrPS 
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where V > is the interaction strength. 

In BCS theory it is assumed that electron pairs have 
zero center-of-mass (angular) momentum, i.e., the pairs 
are condensed in a macroscopic quantum state with 
q = 0, similar to a Bose- Einstein condensate of bosonic 
particles. In the case of a flux threaded ring, Byers and 
YangP, Brenigp-D, and Onsager^ showed that q is gener- 
ally finite and has to be chosen to minimize the kinetic 
energy of the Cooper pairs in the presence of a magnetic 
flux. Nevertheless it is still assumed that pairing occurs 
only for one specific angular momentum q. For conven- 
tional superconductors, this assumption is generally true, 
although for superconductors with gap nodes, the situ- 
ation may be different, as was shown for d-wave pairing 
symmetry in reference^. In this section, we use the as- 
sumption of condensation into a state with one selected 
angular momentum q for all pairs, which allows us to 



write % in the decoupled form 

n = n Q + ^2 [ A * 9 0) c -fc+ g ;c fc1 . + A q {ip)cl t cl k+ql 



V ' 



(7) 



where the order parameter is defined as Aqhp) — 
(V/2) J2k( c kt c ki)- The mean-field Hamiltonian (f7|) is di- 
agonalized with the standard Bogoliubov transformation 



c fcT = u(k)a k+ + v{k)a\_, 
c-k+qi = u{k)a\_ - v(k)a k+ 
with the coherence factors 



u 2 {k) 
v 2 (k) 



1 



E(k,ip) 
e(k,ip) 



(8) 
(9) 

(10) 
(11) 



E(k,<p) 

which depend on ip and q through E(k, ip) = 
^A 2 q + e 2 (k,ip) and e(k, <p) = [e k {<p) + e- k+q (tp))/2. The 
energy spectrum splits into the two branches 

E±(k,(p) = H - 



±JA* + f ?(k,<p), (12) 
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FIG. 7: Eigenenergies E±(k,tp) (121 as a function of flux <p for N = 50 and a fixed order parameter: (a) "large gap" regime 
with A = 0.2 1, (b) "small gap" regime with A = 0.05 1. Blue lines: occupied states, grey lines: unoccupied states. The bold 
line marks the highest occupied state for all ip. In the blue shaded areas in (a), the energy gap has closed due to the Doppler 
shift. 



where the Doppler shift term arises from the different 
energies of the two paired states with momenta k and 
— k + q. The order parameter A q (tp) is determined self- 
consistently from 

i v /(iMfc^))-/(£ + (fc^)) = j fl3l 

2y/A q (<p)* + e*(k,<p) V 

For the discussion of the periodicity of this system, we 
first disregard the self-consistency condition for the order 
parameter and set A q (ip) = A to be constant. Impor- 
tantly, while Hq is strictly <£>o periodic, % is not periodic 
in p if A > 0. The question of periodicity is therefore: 
which periodicity is restored by minimizing E(ip) = {%) 
with respect to q and how is this achieved? 

Figure [6] shows E(ip) for two different values of A for 
q = and q = 1. For small A, E(ip) is still a series of 
parabolae with minima at integer values of <p 7 but the 
degeneracy of the minima is lifted [figure [6] (a) and (b)]. 
For even q, the lowest energy minimum is the one at 
ip = q/2, whereas for odd q, one new minimum at ip = q/2 
emerges. If A exceeds a certain threshold A c , this new 
odd q minimum becomes deeper than the neighboring 
ones [figure [6] (d)]. We have thus identified the second 
class of states with minima in E(ip) at half-integer flux 
values anticipated above and we find that the even and 
odd q minima become equal if A becomes large compared 
to A c , a ring size dependent value which we will deter- 
mine below. It is to be understood that the energies E(ip) 
in figure [6] are not periodic in (p because the q- values are 
fixed, either to q — in (a, c) or to q = 1 in (b, d). In 
loops thicker than the penetration depth, screening cur- 
rents drive the system always into an energy minimum. 
In this case, the flux is then quantized in units of $0/2- 

Let us for the moment assume that the flux value, at 
which the energy minimizing q changes from one integer 
to the next, is well approximated by the half way between 
two minima: q = floor(2(,c + 1/2) (floor(x) is the largest 
integer smaller than x, e.g., p = — > q = 0; p = 1/4 — > 



5 = 1; ip — 3/A^q = 2). Small deviations from these 
values will be discussed for the self-consistent solution in 
section |II C[ The energy spectrum is then $o periodic, 
but discontinuous at the flux values where q changes, 
as shown in figure [7J Clearly distinguishable are now 
the "small gap" (a) and the "large gap" (b) regime: A c 
represents the maximum of the flux-induced shift of the 
energy levels close to Ep, before q changes. If A < A c , 
the energy gap closes at certain values of ip, whereas if 
A > A c , an energy gap persists for all p. 

Although the spectra are <f>o periodic both in fig- 
ure^ (a) and (b), the closing of the gap in the "small gap" 
regime has significant effects on the periodicity of physi- 
cal quantities like E(ip). Even more prominent is the pe- 
riodicity crossover for the persistent current in the ring. 
The supercurrent is given by J (if) — J+(ifi) + J-i}P) — 
(e/h)dE(<p)/d<p, where 

k 

with n + (k) — u 2 (k)f(E + (k,ip)) and — 
v 2 (k)f(E^(k,ipj). J+(p) and J-{<p), as well as J(<p), are 
plotted in figure [8] The contribution J-{<p) forms a $o/2 
periodic saw-tooth pattern, both in the normal and in the 
superconducting state. The $ periodic part in the nor- 
mal state is contained exclusively in J + (p). A flux win- 
dow where E + (k, p) is partially occupied exists in each q 
sector when the energy gap has closed [shaded blue areas 
in figure [7] (a)]. These windows decrease for increasing 
A until J+{p) vanishes for A = A c . In the "large gap" 
regime, the supercurrent is carried entirely by J-(p>) and 
is therefore $o/2 periodic and essentially independent of 
A. The discontinuities in J_ (ip) are not caused by energy 
levels crossing Ep, but by the reconstruction of the con- 
densate when the pair momentum q changes to the next 
integer at the flux values p — (2n — l)/4. Figure [9] shows 
the periodicity crossover of the persistent current in four 
different steps from the $o periodic normal current to 
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FIG. 8: (a) The <J?o periodic persistent current J(</>) (in units of t/$o) m the normal state, (b) The contribution J-(tp) (dark 
blue) is $o/2 periodic and identical in the normal and the superconducting state. The $o periodicity in the "small gap" regime 
is entirely due to J-(tp) (light blue), which vanishes in the "large gap" regime. 
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FIG. 9: Crossover from the $o-P er iodic normal persistent current to the $o/2-periodic supercurrent in a ring with N — 26 at 
T — 0. J(tp) is in units of t / <3?o ■ For this ring size A c ~ 0.24 1. The discontinuities occur where the (^-derivative of the highest 
occupied state energy changes sign. From left to right: A = 0, 0.08 1, 0.16 1, 0.24 1. 



the $o/2 periodic supercurrent. 

Further insight into the current periodicity is obtained 
by analyzing A c . Close to Ep, the maximum energy shift 
is £d = at/2R, and the condition for a direct energy 
gap (or E + (k,p) > for all k, <p) and a $o/2-periodic 
current pattern is therefore A > A c = ep,. The corre- 
sponding critical ring radius is R c = at/2A. It is in- 
structive to compare R c with the BCS coherence length 
£o = hvp/nA, where vp is the Fermi velocity and A the 
BCS order parameter at T = 0. On the lattice we iden- 
tify vp = hkp/m with kp = ir/2a and m — h 2 /2a 2 t, and 
obtain £o = at/ A and thus 2R C = £o- This signifies that 
the current response of a superconducting ring with a di- 
ameter smaller than the coherence length, is generally <£> 
periodical. In these rings the Cooper-pair wavefunction 
is delocalized around the ring. 

We have hereby identified the basic mechanism under- 
lying the crossover from $o periodicity in the normal 
state to < &o/2 periodicity in the superconducting state. 
It is the crossing of Ep of energy levels as a function of 
the magnetic flux that leads to kinks in the energy and 
to discontinuities in the supercurrent (or the persistent 
current in the normal state). If the superconducting gap 
is large enough to prevent all energy levels from cross- 



ing the Fermi energy, the kinks and jumps occur only 
where the pair momentum q of the groundstate changes. 
The latter is true, if the ring diameter is larger than the 
coherence length £ of the superconductor. 

To conclude the discussion of the supercurrent we men- 
tion an issue raised by Little and Parkf^. A simple theo- 
retical model to predict the amplitude of the oscillations 
of T c is the following: For all non-integer or non-half- 
integer values of (p, there is a persistent current J(ip) 
circulating in the cylinder. The kinetic energy E^i n (ip) 
associated with this current is proportional to J 2 (ip), as 
is the energy E(<p) in figure [5] (a). It is therefore sugges- 
tive to subtract -Ekin(</?) from the condensation energy of 
the superconducting state and deduce the oscillations of 
T c from those of EUp). This was done in a first approach 
by Little and Parka^, by TinkharrfED, and by Douglass^ 
within a Ginzburg-Landau ansatz, yet it was later shown 
to be incorrect by Parks and Little in a subsequent ar- 
ticled They wrote that "the microscopic theory [i.e., 
the BCS theory] shows that it is not the kinetic energy 
of the pairs which raises the free energy of the super- 
conducting phase but rather it is due to the differ- 
ence in the energy of the two members of the pairs", i.e., 
e k{p) — e -k+q{ l p)- It is remarkable that the results of 



9 




FIG. 10: Flux threaded annulus with inner radius Ri and 
outer radius 7?2- For a magnetic flux threading the interior 
of the annulus, the radial part of the Bogoliubov - de Gennes 
equations is solved numerically with a discretized radial co- 
ordinate. 



Tinkham and Douglass are nevertheless identical to the 
microscopic resullpS. The notion whether it is the kinetic 
energy of the screening current that causes the oscilla- 
tions, or rather an internal cost in condensation energy 
in the presence of a finite flux, is important insofar as it 
provides an explanation for an intriguing problem: In the 
same way as the pairing of electrons leads to a reduction 
of the fundamental flux period $0 to $o/2, the pairing of 
pairs to quartets would lead to the quarter-period $0/4- 
Then the saw-tooth pattern of the supercurrent becomes 
$0 /4 periodic and the maximum current is only half the 
value for unpaired Cooper pairs. If the oscillation in E(ip) 
were due to the kinetic energy of the pairs, then the for- 
mation of quartets and the <£>o/4 periodicity would be 
energetically favorable. The fact that it is $o/2 periodic 
illustrates the remark by Parks and Little. 



C. Multi channels and self consistency 

Although the one-dimensional ring discussed above 
comprises all the qualitative features of the periodicity 
crossover at T = upon entering the superconducting 
state, some additional issues need to be considered. First, 
the spectrum of a one dimensional ring is special insofar 
as only two energy levels exist that cross the Fermi en- 
ergy in one flux period. This situation is ideal to inves- 
tigate persistent currents, because there are maximally 
two jumps in one period. In an extended loop all ra- 
dial channels contribute at the Fermi energy and have to 
be taken into account. Second, the self-consistency con- 
dition of the superconducting order parameter leads to 
corrections of the results obtained above, and third, the 
periodicity crossover upon entering the superconducting 
state by cooling through the transition temperature T c 
is somewhat different from the T = crossover. These 
points are addressed in this section. 



Here we extend the ring to an annulus with an inner 



radius R\ and an outer radius R2, as shown in figure 10 



For such an annulus, we choose a continuum approach on 
the basis of the Bogoliubov - de Gennes (BdG) equations, 
for which no complications arise from the parity effect. 
For spin singlet pairing the BdG equations are^ 



#nit„(r) = 



E n v a (r) 



2m 



A(r) 



« n (r) 



+ A(r) w„(r), 



2m 



«n(r) 



(15) 



(16) 



+ A*(r) Mn (r), 



with the self-consistency condition (gap equation) for the 
order parameter A(r): 



A(r) = y^ Un (r)<(r)tanh('^ 



(17) 



where V is the local pairing potential. For an annulus 
of finite width we separate the angular part of the quasi- 
particle wavefunctions u n (r), D n (r) using polar coordi- 
nates r = (r, 9) with the ansatz 



(18) 



where k and q are either both even or both odd integers. 
Thus hk is the angular momentum as for the one dimen- 
sional ring and n = (fc, p) with a radial quantum number 
p. The order parameter factorizes into A(r, 9) — A(r)e* 9 
where the radial component 



A(0 = ^ 2^u n (r)<(r)tanh 



2k B T 



(19) 



is real. For a magnetic flux ip threading the interior of 
the annulus we choose the vector potential A(r, 9) = 
eg tp/(2irr), where eg is the azimuthal unit vector. With 



in 

-iV ± -ee 

r 



1 



d r (rd r ) 



1 



(-ido ± V f (20) 



the BdG equations therefore reduce to radial differential 
equations for u n (r) and v n (r): 



E n u n (r) =- 
E n v n (r) 



h 2 d h 2 l 2 
im r 2mr z 



h 2 d r 
2m r 



(rd r )- 



2mr 2 



-p, u n (r)+ A(r)v a (r), 

v n (r)+ A(r)w n (r), 
(21) 



with the canonical angular momenta 

h 

2 
h 

2 



111 u 

III. 



{k + q-2<p), 
(k-q + 2tp). 



(22) 
(23) 
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FIG. 11: Self-consistent calculations for a discretized annulus 
with an inner radius Ri = lOOa^ and an outer radius R% = 
150aj_. The (super-) current (a) J(tp) jumps whenever an 
energy level crosses Ep. The energy levels e(tp) of the normal 
state are indicated by the grey lines, (b) displays the self- 
consistent order parameter A as a function of tp. The lines 
correspond to the pairing interaction V = (orange), V = 
0.28 1 (green), V = 0.32 1 (turquoise), and V = 0.38 1 (blue). 
The black arrows mark the positions of the (/-jump for V = 
0.38t. Here the energy units are t = h 2 /2m e a±. 



For integer and half-integer flux values, equations (21) 
can be solved analytically whereas for an arbitrary mag- 
netic flux a numerical solution is required (for details, 
see reference^). Within this procedure, the radial coor- 
dinate is discretized into N± values r n separated by the 
distance aj_ = (i?2 — Ri)/N±. The number q plays the 
same role as in the previous section. Here we choose q for 
each value of the flux to minimize the total energy of the 
system. The flux for which q changes to the next integer 
can therefore deviate from the values (2n— 1)/4, at which 
we fixed the jump to the next q for the one-dimensional 
model. 

In the normal state (A = 0), the number of eigenstates 



sufficiently close to Ep which may cross Ep as a function 
of tp is controlled by the average charge density n and is 
approximately (R2 — Ri)/a± for n = 1. For each cross- 
ing, a jump appears in the current as a function of tp, as 
shown in figure [IT] (a) . The persistent current is there- 
fore proportional to the level spacing at E-p for each flux 
value, i.e., it is small for a large density of states and 
large for a small density of states and vanishes in the 
limit of a continuous density of states (c.f. reference^). 
The amplitude of the normal persistent current is thus 
a measure for the difference of the energy spectrum at 
integer and half integer flux values. It is maximal for 
the one dimensional case, but might be very small in real 
metal loops (c.f. measurements of the Aharonov-Bohm 
effect in metal ring 

Upon entering the superconducting state, an energy 
gap develops around Ep preventing energy levels from 
crossing Ep. Thus the persistent supercurrent arises as 
in the one dimensional ring independent of the density 
of states. The large jumps in the supercurrent appear at 
the value of tp where the energies of the even-g and odd- 
q states become degenerate and q switches to the next 
integer. In the flux regimes with no crossings of Ep the 
supercurrent is linear and the total energy quadratic in 
ip. For the largest value shown (A = 0.006/;), there is a 
direct gap for all values of tp. Even for this large A, the 
current and the energy are not precisely $o/2-periodic 
because of the energy differ ence of the even and odd q 
states in finite systemJ^^H. The offset of the g-jump is 
only relevant for values of the pairing interaction V for 
which A is finite for all tp. In figure 11 (b), the offset 
is clearly visible for the largest two values of V (marked 
with black arrows). Its sign depends on the shape of the 
annulus and the value of V — the offset changes sign for 
increasing V (cf. reference^). 

The introduction of self-consistency for the order pa- 
rameter does not fundamentally change these basic obser- 
vations [figure [TTj (b)]. One finds that A(i?i) < A(i? 2 ), 
but if (i?2 — R1J/R1 % 1, the difference is small. In the 
following, we denote the average of A(r) by A. The 
crossover is then controlled by the pairing interaction 
strength V, for which we chose such values as to repro- 
duce the crossover from the normal state to a state with 
direct energy gap for all flux values. The order parame- 
ter A is now also a function of tp. If A(tp = 0) < 0.006 1 
[c.f. figure 11 (b)], the gap closes with tp, and A decreases 
whenever a state crosses Ep. Unlike in a one dimensional 
ring, A does not drop to zero at the closing of the en- 
ergy gap, but decreases stepwise. This is because in two 
or three dimensions, A is stabilized beyond the depair- 
ing velocity by contributions to the condensation energy 
from pairs with relative momenta perpendicular to the 
direction of the current flow; the closing of t he ind irect 
energy gap does not destroy superconductivity^^^. 

Experimentally more relevant is to control the 
crossover through temperature. With the pairing in- 
teraction V sufficiently strong to produce a T = en- 
ergy gap much larger than the maximum Doppler shift, 
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FIG. 12: The order parameter A(tp) and the persistent current J(<p) for the temperature driven transition from the normal to 
the superconducting state in an annulus with inner radius Ri — 30a± and outer radius R2 = 36a±. The pairing interaction is 
Vo = 0.7t, with a critical temperature of ksT c « 0.0523t for zero flux. For these parameters A(T = 0) « O.lt. The lines (from 
top to bottom) correspond to the temperatures feT = 0.0513i (blue), feT = 0.0520i (turquoise), and feT = 0.0522t (green). 
Notice that A is slightly different for the flux values ip = and tp — ±1/2. 



the crossover regime is reached for temperatures slightly 
below T c . For the annulus of figure [Toj the crossover 
proceeds within approximately one percent of T c . The 
crossover regime becomes narrower for larger rings, pro- 
portional to the decrease of the Doppler shift. In the 
limit of a quasi one-dimensional ring of radius R we can 
be more precise: If we define the crossover temperature 
T* by A(T*) = A c and assuming A c -C A, we can use 
the Ginzburg-Landau form of the order parameter 



A(T) 
A(0) 



1.75W1 - 



T 



(24) 



and obtain 

T r — T* 



3.1A 2 (0) 12.4A 2 (0)(i?/a) 2 



El 



3.l(k B T c y(R/af 



(25) 



where we used the relations A c = at/2R and Ep = 2t 
for a discretized one-dimensional ring at half filling, and 
the BCS relation A 2 (0) « 3.1(fc B T c ) 2 . For a ring with a 
radius of 2500 lattice constants (« 10 |xm) and A(0) = 
0.01 1 (w 3meV) one finds the ratio (T c - T*)/T c » 
1.3 x 10~ 4 . This is in reasonable quantitative agre ement 
with the experimental results of Little and Parks^^, 
discussed also by Tinkham in reference^. Their pre- 
diction was similar to equation ( 25 ) , up to a factor in 



which they include a finite mean free path. But they did 
not included the difference introduced through even and 
odd q states. This difference was considered in calcula- 
tions of T c by Bogachek et aZP^in the single-channel limit 
and found to be exponentially small. A detailed study 
of the normal- to superconducting phase boundary was 
also done by Wei and Goldbart in reference^ in which 
they considered the $0 periodic contributions. In equa- 



odd q. Although quantitative predictions of T c — T* of 
the theory presented here might be too large as compared 
to the experiment, it serves as an upper limit because it 
describes the maximum possible persistent current. In- 
homogeneities and scattering processes in real systems 
further reduce the difference of the energy spectra in the 
even- and odd-g flux sectors and thereby reduce T c — T* . 

For temperatures close to T c , the difference of the 
eigenenergies of even and odd q states is less important 
than at T = 0. Thus the deviation from the $o/2 period- 
icity of the current and of the order parameter is smaller. 
Furthermore, persistent currents in the normal state are 
exponentially small compared to the supercurrents below 
T c . Their respective $0 periodic behavior is therefore es- 
sentially invisible for the flux values where A = 0. In fig- 
ure 12 the difference between A(ip = 0) and A(ip = 1/2) 



tion ( 25 1 the value of A(0) is in fact different for even and 



is still visible, but the corresponding differences in the 
current are too small. Only for a superconductor with 
very small T c , we expect the periodicity crossover to be 
visible. 

Although we found within the framework of the BCS 
theory, that the crossover to a $o/2 periodic supercurrent 
takes place slightly below T c , detailed studies by Ambc- 
gaokar and EckerrP^ and by von Oppen and RiedeP^ 
including superconducting fluctuations reveals that the 
crossover might actually take place above T c . This fluctu- 
ation driven crossover is broader than the BCS crossover 
with a similarly, exponentially suppressed $0 periodic 
normal current contribution. For a superconductor with 
a T c small enough to observe a normal persistent cur- 
rent above T c , Eckcrn and Schwab suggested that the 
crossover regime, where both $0 and $o/2 periodic cur- 
rent contributions ar e pres ent, should be observable at a 
temperature T w 2Tp^. 

The discussion of the periodicity crossover in a multi- 
channel loop also gives insight into the flux periodicity of 
loops of unconventional superconductors with gap nodes 
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like a d-wave superconductor. In nodal superconduc- 
tors the density of states is finite arbitrarily close to Ep. 
Therefore some energy levels cross Ep as a function of 
the flux, regardless of the size of the order parameter and 
consequently, the "small gap" situation extends to arbi- 
trarily large loopsP^l Of course, the number of energy 
levels crossing Ep decreases with increasing ring size and 
thus also the $0 periodic contribution to the supercur- 
rent. The dependence of this <&o periodic contribution on 
the ring size depends on the order parameter symmetry. 
The careful study in reference^ revealed that for d-wave 
superconductors, the relation between the $0 and $o/2 
periodic current contributions is proportional to l/Ri- It 
was estimated that for a ring of a cuprate superconduc- 
tor with a circumference of ~ 1 u\m, this ratio is about 
1% and should be observable experimentally. 



III. CONCLUSIONS 

We analyzed the crossover from the $0 periodic persis- 
tent currents as a function of magnetic flux in a metallic 
loop to the $0/2 periodic supercurrent in the groundstate 
of the loop. We considered conventional s-wave pairing 
in a one-dimensional as well as in a multi-channel annu- 
lus. Although a one-dimensional superconducting ring 
is a rather idealized system, it proves valuable for dis- 
cussing the physics of this crossover, which includes the 
emergence of a new minimum in the free energy for odd 
center-of-mass angular momenta hq of the Cooper pairs 
and the restoration of the flux periodicity of the free en- 
ergy. The physical concepts, which we illustrated in a 



simplified form in section II B remain thereby valid even 
in the more complex context of the self consistent calcu- 
lations on the annulus. 

In the superconducting state, a distinguished minimum 
in the free energy develops at ip — q/2. Choosing the 
proper value for q at each flux value leads to a series 
of minima at integer and half-integer flux values which, 
however, differ in energy for finite systems. In rings with 
a radius smaller than half the superconducting coherence 
length, the two electrons forming a Cooper pair are not 
forced to circulate the ring as a pair, and the supercurrent 
shows a $0 periodicity. Only if the order parameter A is 
larger than the maximal Doppler shift ep,, the supercur- 
rent is $o/2 periodic. This is equivalent to the condition 
that the maximum flux induced current is smaller than 
the critical current J c . Assuming that the relations ob- 
tained from the one-dimensional model remain valid on 
a ring with finite thickness R2 — R\ <C R\ , as indeed sug- 
gested by the multi-channel model, the critical radius to 
observe $o/2 periodicity, R c = at/2A, would be of the 
order of 1 |i.m for aluminum rings. Within the temper- 
ature controlled crossover upon cooling through T c , $ 
periodicity might by difficult to observe since the differ- 
ences in the energy spectra for integer and half-integer 
flux values are exponentially suppressed by temperature. 
$0 periodicity is therefore only observable if a normal 



persistent current would be observable at the same tem- 
perature if superconductivity was absent. 

In the introduction we referred to experiments where 
flux oscillations with "fractional periodicities", i.e., frac- 
tions of $o/2, were observed. Among various suggested 
origins, there is one particularly elegant approach based 
on a standard two-electron interaction. Consider the or- 
der parameter A q (ip) for electron pairs with center-of- 
mass angular momentum hq. In real space, q describes 
the phase winding of the order parameter A(9, ip) = 
A(ip)e lq6 , where 9 is the angular coordinate in the ring 
and A(ip) is real. To ensure that A(9, ip) is a single valued 
and continuous function, q must be an integer number. 
If, however, A(9, ip) is zero somewhere on the ring, it can 
change sign. Such a sign changing order parameter is 
modeled as 



A(d, l p)= 1 -[A ( l p) + A q (tpy« e ] 



^AMe^/ 2 cos(|0), (26) 



which displays a phase- winding number q/2 if A q — A , 
and consequently a vanishing supercurrent at the frac- 
tional flux value g$o/4. In momentum space, A(9,ip) 
is represented by the two-component order parameter 
{A ((p), A q (<p)}. Such a superconducting state is typi- 
cally referred to as a "pair-density wave" (PDW) state^, 
since the real-space order parameter is periodically mod- 
ulated and therefore q can no longer be interpreted as an 
angular momentum. Agterberg and Tsunetsugu showed 
within a Ginzburg-Landau approach, that a PDW super- 
conduct or ind eed allows for vortices carrying a $o/4 flux 
quantuirPE] 

Based on a microscopic model, it was shown in ref- 
erence^ that a PDW state cannot result from an on- 
site pairing interaction. However, the PDW state can be 
a stable alternative for unconventional superconductors 
with gap nodes, specifically for d-wave superconductivity 
as realized in the high-T c cuprate superconductors. More 
work is, however, needed to analyze under which condi- 
tions the PDW state indeed develops an energy minimum 
at multiples $0/4, and these minima become degenerate 
in the limit of large loops. The notion of the absence of 
coexistence of Cooper pairs with different center-of-mass 
momenta for conventional superconductors (as, e.g., in 
the PDW state) justifies the reduction of the sum over 
q in the Hamiltonian Q in section II B to one specific q 
in order to derive the periodicity crossover in supercon- 
ducting rings. 

Here we mention a further system where a similar 
mechanism as described above may lead to $0 /4 flux pe- 
riodicity: Sr2Ru04. Experimental evidence exists that 
Sr 2 Ru04 is a spin triplet p-wave superconductor. A 
triplet superconductor can be represented by the two- 
component order parameter {A]^(ip), A^(ip)}, with the 
center-of-mass momenta q\ and qi for the s z = 1 and the 
s z = —1 condensates. This realizes a similar situation 
as for the PDW state where $o/4 periodicity is possi- 
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